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Abstract 

We compute the contribution of kinks on cosmic string loops to stochastic 
background of gravitational waves (SBGW). We find that kinks contribute at 
the same order as cusps to the SBGW. We discuss the accessibility of the total 
background due to kinks as well as cusps to current and planned gravitational 
wave detectors, as well as to the big bang nucleosynthesis (BBN), the cosmic 
microwave background (CMB), and pulsar timing constraints. As in the case 
of cusps, we find that current data from interferometric gravitational wave 
detectors, such as LIGO, are sensitive to areas of parameter space of cosmic 
string models complementary to those accessible to pulsar, BBN, and CMB 
bounds. 



1 Introduction 



Topological defects are remnants of spontaneously broken local or global sym- 
metries. The simplest and the most well-known example of the former one 
is the Abrikosov-Nielsen-Olesen flux tube [1], which originates from sponta- 
neously broken U(l) gauge symmetry. Most of the attention in the literature 
has been focused on defects originating from broken gauge symmetries, since 
grand unified theories have gauge symmetries which are eventually sponta- 
neously broken down to the symmetry of the Standard Model. Cosmic strings 
are one dimensional topological defects predicted by a large class of unified 
theories [2, 3, 4]. Cosmic strings were first considered as the seeds of struc- 
ture formation [5, 6], however, later, it was discovered that cosmic strings 
were incompatible with the cosmic microwave background (CMB) angular 
power spectrum. Cosmic strings can still contribute to structure formation, 
but they cannot be the dominant source. Cosmic strings are also candidates 
for the generation of other observable astrophysical phenomena such as high 
energy cosmic rays, gamma ray burst and gravitational waves [3, 7, 8, 9]. 
Furthermore, recently it has been shown that in string-theory-inspired cos- 
mological scenarios cosmic strings may also be generated [10]. They are 
referred to as cosmic superstrings. This realization has revitalized interest in 
cosmic strings and their potential observational signatures. There are some 
important differences between cosmic strings and cosmic superstrings. The 
reconnection probability is unity for cosmic strings [3, 11]. Cosmic super- 
strings, on the other hand, have reconnection probability less than unity. 
This is a result of the probabilistic nature of their interaction and also the 
fact that it is less probable for strings to meet since they can live in higher 
dimensions [12]. The value of p ranges from 10~ 3 to 1 in different theories 
[13]. Cosmic superstrings could also be unstable, decaying long before the 
present time. In this case, however, they may also leave behind a detectable 
gravitational wave signature [14]. 

In the early universe, a network of cosmic strings evolves toward to an 
attractor solution called the "scaling regime". In the scaling regime the 
statistical properties of the network, such as the average distance between 
strings and the size of loops at formation, scale with the cosmic time. In 
addition, the energy density of the network remains a small constant frac- 
tion of the energy density of the universe. For cosmic superstrings in the 
scaling regime, the density of the network p is inversely proportional to the 
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reconnection probability p, that is p oc p ? . The value of (3 is still under 
debate [15, 16, 17], and as a placeholder in our analysis we assume that 

= 1- 

The gravitational interaction of strings is characterized by their tension 
fi, or more conveniently by the dimensionless parameter Gpi, where G is New- 
ton's constant. The current CMB bound on the tension is Gfi < 6.1 x 10 -7 
[18, 19]. It was first believed that gravitational radiation from cosmic strings 
with Gfj, <C 10 7 would be too weak to observe. However it was later shown 
that gravitational radiation produced at cusps, which have large Lorentz 
boosts, could lead to a detectable signal [20, 21, 22]. Gravitational radiation 
bursts from (super) strings could be observable by current and planned grav- 
itational wave detectors for values of G\i as low as 1CT 13 , which may provide 
a test for a certain class of string theories [23]. Indeed, searches for burst 
signals using ground-based detectors are already underway [24]. 

A gravitational background produced by the incoherent superposition of 
cusp bursts from a network of cosmic strings and superstrings was considered 
in [25]. In this paper we extend this computation to include kinks, long- 
lived sharp edges on strings that result from intercommutations, and find 
that kinks contribute at almost the same level as cusps. We investigate the 
detectability of the total background produced by cusps and kinks by a wide 
range of current and planned experiments. A similar calculation for the case 
of infinite strings has been undertaken in the recent paper [26], see also [27]. 

The organization of the paper is as follows: In Sect. 2 we consider gravi- 
tational waves generated by cusps and kinks in the weak field limit [28]. In 
this section we follow the conventions of [20, 21], and more details can be 
found in these references. In Sect. 3 we derive the expression for the stochas- 
tic background, which is a double integral over redshift and loop length. In 
Sect. 4 we evaluate integral analytically with certain approximations, which 
results in a flat distribution for larger values of the frequency. Finally in Sec. 
5 we numerically evaluate the background and discuss the observability by 
various experiment. 
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2 Gravitational Radiation 



In this section we consider gravitational waves created by cusps and kinks. 
For completeness we follow closely the analysis in [20, 21], and reproduce a 
number of their results. We begin with a derivation for the metric pertu- 
bation in terms of the Fourier transform of the stress energy tensor of the 
source. We then write the stress energy tensor for a relativistic string and 
compute its Fourier transform. Using these results we then compute the 
gravitational waveforms produced by cusps and kinks on cosmic strings. 

2.1 Calculation of metric perturbations 

Gravitational waves from a source can be calculated using the weak field 
approximation [28], 

9iiv = Vt^ + h^u, (1) 

where r]^ is the Minkowski metric with positive signature and is the 
metric perturbation. In the harmonic gauge, g^ u T x u = 0, the linearized 
Ricci tensor is 

Rp* ^ \d x d x h, K . (2) 
Substituting into Einstein's equations yields 

Rfiu - ^g^R ^ ^(<9 A <9 A /v - ^V^d x d x h) = -8nG% u , (3) 

where R is the Ricci scalar, T^ u is the energy momentum tensor of matter 
and h = r\ tlv U lv . Defining h^ v = — ^rj^h further simplifies Eq. (3), 

d x d% u = -16nG% u , (4) 

which is a wave equation with a source term. We can rewrite this equation 
in the frequency domain as, 

(w 2 + V 2 ) V(f, w) = -16nG% u (x, w), (5) 

where 

h^(x : w) = J dte iw %„(x,t). (6) 
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Eq. (5) can be solved by using the Green's function for the operator 
w 2 + V 2 , which is 

giro | x— x' | 

g(x-x,w) = — — — -. (7) 
\x — x | 

Therefore metric perturbations are given by 

h^Xjiv) = — 16nG J d 3 x'Q(x — x', w)T^(x', w) 

= -16nG^-—% u {k,w), (8) 



\x\ 



where k = wx and 



1 T 

%„(k,w) = -jf dtJd*Jj^- % ^%,{?,t), (9) 

where T is the fundamental period of the source. Eq. (8) relates energy 
momentum tensor to gravitational waves. The next step is to calculate the 
energy momentum tensor of cusps and kinks on cosmic strings. 

2.2 Energy Momentum Tensor of Cosmic Strings 



In the thin wire approximation, the dynamics of strings is described by the 
Nambu-Goto action [3, 7] 

S = -fiJ drdo^i, (10) 

where a and r are world-sheet coordinates and ji is the string tension. 7 is 
the determinant of the induced metric 

lab = V^9 a X' i d b X v , (11) 

where a and b denote world sheet coordinates. The equation of motion 
following from Eq. (10) is 

(d 2 T - d 2 a )X» = 0, (12) 
The solution must also satisfy Virasoro conditions 

X ■ X + X' ■ X' = and X • X' = 0, (13) 



4 



where dot and prime denote derivatives with respect to r and a 
respectively. If we define o± — r ± a, the equation of motion becomes 

+ d_X» = 0, (14) 

which is solved by left and right moving waves, 

X» = \{X»{a + ) + Xt{*-)) . (15) 

Furthermore Virasoro conditions in Eq. (13) simplify to 

±1 = 1, (16) 

where dot now represents the derivative with respect to the (unique) 
argument of the functions X±. We require that X^(cr, r) is periodic in a 
with period /, which is the length of the loop. This implies that the 
functions X± are periodic functions with the same period. The period in t 
is 1/2 since X»(a + 1/2, r + 1/2) = X»(a, r). 

The energy momentum tensor corresponding to the Nambu-Goto action can 
be calculated by varying Eq. (10) with respect to the metric, which yields 

5S 



: -2-p- = n [ drda(X^X v - X' U X^) 5 (4) (a; - X) 
= | / do-do+{XlX v _ + X v _Xt) S^(x - X). (17) 

Inserting this expansion into Eq. (9) gives us the energy momentum tensor 
in momentum space 

%M = ^ f da.da + X^X%-i^ k - x -\ (18) 



where we define 

X^X V J = ^(X»X»+X»X»). (19) 

The nice property of Eq. (18) is that two integrals can be calculated 
independently, 

l£(k) = f da ± X^e~^ x± , (20) 
Jo 
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and the energy momentum tensor can be expressed in terms of I± as 
follows; 

%„(k) = (21) 

where we used TJ = |. In the following subsection we calculate 7± for cusps 
and kinks. 

2.2.1 Cusps 

Let us start with the geometrical interpretation of Eq. (16). It tells us that 
X± trace a unit sphere centered at the origin, which is called Kibble- Turok 
sphere. Integrating X± and using the periodicity, we get 

[ l X ± (a±)da ± = 0, (22) 
Jo 

which implies that X± cannot lie completely in a single hemisphere and 
therefore they intersect at some point (s). We choose our parametrization 
and the coordinate system such that the intersection occurs at the 
parameters a± — at the origin, that is X±(0) = 0. X±(a±) and X±(a±) 
can be expanded around a± — 

X£(a ± ) = l^ ± + l -X^l+ l -X { i^al (23) 
X£(a ± ) = F ± + X£a ± + ±X^al. (24) 

where l± = X±(0). We can easily find the shape of X± at r = (a± = ±cr), 

X"(v,r = 0) = 1{X$(&) + Xtira)) 

= -(X»+X»)a 2 + —(X^ + X { ^)a 3 . (25) 

In order to visualize the shape of the string around the origin, we can 
choose the coordinate system such that (X + + X_) lies on the and 
define x = \\X+ + X_|cr 2 . Let us also denote the direction of X^ + X^ by 
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y, which is not necessarily orthogonal to x. If we define 

y = -^[X^^ + X^^\a 3 , we see that y oc x%, which has a sharp turn at 

x = 0, which is referred to as cusp. 

We can calculate I± for cusps using the expansion in Eq. (23). First of all, 
we note that the first term Eq. (24) is pure gauge, it can be removed by a 
coordinate transformation. Furthermore imposing Virasoro condition in Eq. 
(16) gives 

l± ■ X± = 0, and 1± • Xi 3) = -X|. (26) 

When the line of sight k is in the direction of / we have k = wl, which gives 
—i k ■ X± = ^wX±a±. If we plug in the expansion in Eq. (23) into Eq. (20) 
we get, 

im = n fdaae^ 3 = 2 f X * ^ (27) 

J ° 31X1/3) Ol^l) 73 

Replacing w with 2nf gives 

im = C3J/-i, (28) 
%M = ^\f\-H"C- } (29) 



/i - (327r/3) 1 / 3 

|X±|3 

Since X± is periodic with period /, X expanded as 



where C± = i ^ r(i/3) , ,-> TI • Fi na % we nee d to estimate |X±| = |X±| 



±(a ± )=J2c n e^ n ^, (30) 

n 

where the expansion coefficients c n are constrained by |X±| = 1. If the 
string is not too wiggly, c n is nonvanishing for only small n, therefore we 
can estimate \X±\ ~ ^f-. Combining all the pieces together and neglecting 
decimal points in the numerical coefficient, we express the trace of the 
metric perturbations as 

h^f) = H\ = ^\f\-l. (31) 

We can express r as a function of z 

1 f z dz' 1 . . . . 

= —<p r (z), (32) 



H Jo n(z') H 
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where H Q is the Hubble constant today and %{z) is the Hubble function 
given by 

U{z) = (n M {\ + zf + n R {l + zf + Q A ) 1/2 . (33) 

The numerical values for the constants in this equation are Qm — 0.25, 
tt R = 4.6 x 1(T 5 , n A = 1 - tt R - VL M and H = 73km/s/Mpc. 
Note that / in Eq. (31) is the frequency of the radiation in the frame of 
emission. In order to convert it to the frequency we observe today, the 
effect of the cosmological redshift must be included. The frequency in the 
frame of emission, /, is related to the frequency we observe now, / now > by 
the relation / = (1 + z)f now . After redshifting properly 1 , Eq. (31) becomes 

h>^f,zj)= G ^ li \f\-i, (34) 

(1 + Z)a<p r {z) 

where we dropped the subscript now. 
2.2.2 Kinks 



Calculation of kink radiation is similar to the cusp case. The form of 7+ is 
the same as the cusp result. Ift has a discontinuity at the cusp point and 
needs a different treatment. Let us describe the kink (at o"_ = and 
X± = 0) as a jump of the tangent vector from 1% to 1% . At the first order 
one can replace approximate by 1% for er_ < and 1% for cr_ > 0. At this 
approximation, one gets 

I-(k) = [ l/2 da-XHe-i"- x - ~ - ( -?L - , (35) 

where we dropped two oscillatory terms. The exact value of Eq. (35) 
depends on the sharpness of the kink, l\ ■ l 2 [29] , however we will assume 

1 One should note that replacing / in Eq. (31) with (1 + z)f now is not correct, since this 
replacement will scale the argument and the amplitude of (/) by a factor of j4^, which 
is the reflection of the fact that the measure of Fourier integral is not dimensionless. Since 
redshifting should change the argument but not the amplitude, one needs to multiply the 
result by 1 + z so that the amplitude remains the same. Equivalently, one can define 
Logarithmic Fourier Transform, as discussed in Ref. [20] , such that the measure of the 
transform becomes dimensionless. 
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that the average value of this quantity is of order one. Combining this result 
with i+ we get the frequency distribution of the radiation from a kink as 

h™(f,z,l) = - G ^ H \ f^ 3 . (36) 

(1 + Z)3(p r ( Z ) 

It is important to note that in the derivation of Eqs. (34) and (36) we 
assumed that the line of sight is in the direction of the motion of the 
cusp or kink, Z M . It is easy to show that I± (Eq. (20)) decay exponentially 
with the angle between k and 1[21] . Therefore Eqs. (34) and (36) are valid 
for angles smaller than 

e m = r - (37) 

We implement this condition with a O-function in the amplitude. 



3 Stochastic Background 

The stochastic gravitational background [25] is given by 

fW/) = *f>f 4, J (38) 

where h(f, z, I) is given in Eqs. (34) and (36) and d j^'^ is the observable 
burst rate per length per redshift, which will be defined below. We take the 
number of cusps (kinks) to be one per loop. If we define the density (per 
volume) of the loops of length I at time t as n(l,t), the rate of burst (per 
loop length per volume) can be expressed as ^rp, where 1/2 factor is the 
fundamental period of the string. However, this is not the observable burst 
rate since we can observe only the fraction of bursts that is beamed toward 
us. Including this fraction we obtain 

//nVrOll • :) '^A(,/./|. (39) 



dldz u r v /y I 
where (1 + z)^ 1 comes from converting emission rate to observed rate, and 
HQ S {py(z) follows from converting differential volume element to the 
corresponding function of redshift z, 

dV = *Ka\tydr = ^f^f^ dz ee H^v(z)dz, (40) 
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where a(t) is the cosmological scale factor. A(z, f, I) is the fraction of the 
bursts we can observe. Geometrically the radiation from a cusp will be in a 
conic region with half opening angle 9 m (Eq. (37)) and outside the cone it 
will decay exponentially. To simplify the calculation we assume that the 
radiation amplitude vanishes outside this conic region, which will be 
implemented by a O-function. We can express the corresponding solid angle 
in terms of the opening angle by using the following relation 

ft m = 27T(l-COS0 m )~7<. (41) 

Thus the probability that the line of sight is within this solid angle is 

£ - C/4, (42) 

which is referred to as the beaming fraction of the cusp. We combine the 
cutoff for large angles and beaming effect into 

A(*, /, « ^ 9(1 - 6 m (z, f,l)). (43) 

It is important to note that cusps are instantaneous events, and it is 
possible to observe their radiation only if the line of sight happens to be 
inside the cone of radiation. The beaming fraction, Eq. (42), which is 
proportional to d 2 ml is the fraction of the time the line of sight is inside the 
cone of radiation. In contrast, kinks radiate continuously-as kinks travel 
around a string loop they radiate in a fan-like pattern. Therefore radiation 
cone of a kink will sweep a strip of width 26 m and an average length 7r on 
the surface of the unit sphere as it travels around the cosmic string loop. 
That is, the probability of observing radiation from a kink is 

4tt 4tt 2 ' V ; 

For kinks the cutoff for large angles and beaming factor that enters the rate 
is therefore 

A<*>(s, /, « 9m{Z 2 Ll) Q(l - 9 m (z, /,/)). (45) 

Inserting this result into Eq. (38) gives the background radiation Q gw (f) as 
a double integral over / and z, which needs to be evaluated numerically. 



10 



Finally we need to discuss the form of the loop density, n(l,t) in Eq. (39). 
To do this, it is convenient to first convert the cosmic time t to a suitable 
function of redshift z using the following relation 



dz 
~dl 



(l + z)H H(z), 



(46) 



which can be integrated to give 



t = H l / 

J z 



oc 



dz 1 



(47) 



(i + z')n(z') 



Below we discuss the two main contending scenarios for the size of cosmic 
string loops. 

3.1 Small Loops 

Early simulations suggested that the size of loops was dictated by 
gravitational back reaction. In this case the size of the loops is fixed by the 
cosmic time t, and all the loops present at a cosmic time t, are of the same 
size a t. The value of a is set by the gravitational back reaction, that is 
a oc TGfi (In Sect. 5 we parameterize a by a = eYGji where e is a 
parameter we scan over.) The constant V is the ratio of the power radiated 
into gravitational waves by loops to Gfi 2 . Numerical simulation results 
suggest that T ~ 50. Therefore the density is of the form 



where p is the reconnection probability. The overall coefficient is estimated 
by simulations (for a review see [3]) which show that the density in the 
radiation domination era is about 10 times larger the one in the matter 
domination era. This behavior of the density can be implemented by a 
function, c(z), which converges to 10 for z ^> z eq and to 1 for z z eq . 
Therefore the density can be written as 



n(l,t) oc (pFGii)-H- 3 5(l - at), 



(48) 



n(l,t) = c(z)(pTGfi)-H- 3 5(l - at), 



(49) 



where [20] 



c(z) = 1 + 



9z 



(50) 



z + z, 
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Such a distribution simplifies the calculation of SBGW since the /-integral 
in Eq. (38) can be evaluated trivially to yield 

^ 2 r^r»u2,*„H*RM 



*W/) = ^2 Jdzjdlh 2 (f,z,l) 



3Hq J J v "' ' ' dzdl 
2cGpm 2 H 1 /3 r c(z) Vv Q(l-[f(l + z)a Vt ]- 1/3 ) 

dz t~. zzz - 10/3 • y bl ) 



For kinks, we have a similar integral , 

QK(f) _ acg^hT f . A^y® i 1 - [/a + *wr 1/3 ) 

We analytically evaluate the integrals in Eqs. (51) and (52) in Sec. 4 with 
certain approximations, and perform numerical integration in Sec. 5. 



3.2 Large Loops 

Recent simulations [30, 31, 32] suggest that the size of the loops is set by 
the large scale dynamics of the network, and that the gravitational 
back-reaction scale is irrelevant. In Ref. [32] is found that the loop 
production functions have peaks around a ~ 0.1, which is the value we use 
below (for large loop case). For long-lived loops, the distribution can be 
calculated if a scaling process is assumed (see [3]). In the radiation era it is 

n(l,t) = x r H(Z + rG)rf)-2, 

l<at, t<t teq (53) 

where Xr ~ OA^a 1 ^ 2 , and £ is a parameter related to the correlation length 
of the network [22]. The numerical value of £ is found in numerical 
simulations of radiation era evolution to be about 15 (see Table 10.1 in [3]). 
The upper bound on the length arises because no loops are formed with 
sizes larger than at. For t > t eq (the matter era) the distribution has two 
components, loops formed in the matter era and survivors from the 
radiation era. Loops formed in the matter era have lengths distributed 
according to, 

ni(l,t) = X mt- 2 (l + TGi2t)- 2 , 

at teq - TGfi(t - t t J <Kat, t> t Uq (54) 
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with Xm ~ 0.12^, with ( 4 (see Table 10.1 in [3]). The lower bound on 
the length is due to the fact that the smallest loops present in the matter 
era started with a length a t eq when they were formed and their lengths 
have since decreased due to gravitational wave emission. Additionally there 
are loops formed in the radiation era that survive into the matter era. 
Their lengths are distributed according to, 

n 2 (l,t) = Xrtlft-^l + TGfityl, 

l<at teq -TGfi(t-t teq ), t>t teq , (55) 

where the upper bound on the length comes from the fact that the largest 
loops formed in the radiation era had a size a t eg but have since shrunk due 
to gravitational wave emission. 

The cusp spectrum has been calculated in [25] and the result shows that 
the spectrum is flat for larger 2 values of /. Later we will show that this is 
also the case for kink spectrum. This is rather unexpected since £l(f) has 

4 1 

an explicit /~s and /~5 dependence for cusps and kinks, respectively. The 
only other / dependence comes from the 9 functions. In the following 
section we show analytically that the / dependence coming from the 

4 1 

function is of the form /5 and fs for cusps and kinks respectively so that 
the spectrum is indeed flat for large values of the frequency /. 
Before we start calculating the SBGW, we should mention a crucial 
observation due Damour and Vilenkin [20] . SBGW generated by a network 
of cosmic strings includes bursts which occur infrequently, and the 
computation of Q gw (f) should not be biased by including these large rare 
events (i.e. events with low rate). If the loop density is taken of the form 
given in Eq. (48), the rate is specified by the redshift only. Therefore the 
condition on the rate can be implemented by a cutoff on redshifts such that 
large events for which the rate is smaller than the relevant time-scale of the 
experiment are excluded (see Eq. (6.17) of [20]). However, when loops are 
large the situation is more complicated because at any given redshift there 
are loops of many different sizes given in Eqs. (53) and Eq. (54). This case 
has been dealt with in [25] as follows: instead of integrating over the 
variables / and z in Eq. (38) one integrates over h and z where h is defined 

2 In the following section we show that the spectrum is flat for / » H °V^ f or sma n 
loops and for / » ^Gur^ ^ or l ar § e loops. 
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in Eqs. (34) and (36) and imposes the cutoff limit on the h integral. The 
cutoff is defined as 

r dk i dz Pk = /> ( 5e ) 

Jh* J azah 

where = tOtJI- Eq. (56) is solved for h* and used to exclude rare 
event using the following integral (instead of Eq. (38)) 

^'m/r^hzk- (57) 

This procedure removes large amplitude events (those with strain h > h*) 
that occur at a rate smaller than /. Fig. 1 shows the spectrum for kinks 
and cusps for small loops. For the top curves (red and green) we have, 
Gfi = 2 x 10~ 6 , p = 10 -3 and e = 10 -4 , whereas for the bottom two curves 
(blue and pink) G\i = 10~ 7 , p = 5 x 1CT 3 and e = 1 ( e = pg- ). 



Cusp'(1) 




,10 



f (Hertz) 

Figure 1: Kink and Cusp spectrum for small loops: (1) Gfi = 2 x 10 
p = 10" 3 and e = 10" 4 , (2) Gfi = 10^ 7 , p = 5x 10" 3 and e = 1 . 

Fig. 2 shows the spectrum for large loops. For the top curves (blue and 
pink), which are almost identical, we have, Gfi = 10~ 7 and p = 5 x 10 -3 , 
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whereas for the bottom two curves (red and green) Gfi = 10 9 and 

p = 5x 1CT 2 . Here we note that for / ^> the spectrum is flat for both 




f(Hertz) 



Figure 2: Kink and Cusp spectrum for large loops: (1) Gfi = 10 7 and 
p = 5x 10" 3 , (2) Gfi = 10" 9 and p = 5 x HT 2 . 



cusps and kinks. 



4 Analytical Approximation for the 
Stochastic Background 

In this section we evaluate the spectrum analytically and show that the 
spectrum is constant for large values of /. Our main goal is the discuss the 
dependence of the spectrum on the parameters: Gfi, e = and p for 
small loops and G\x and p for large loops. We limit our discussion to large 
values of /, for which the spectrum gets the dominant contribution from 
the loops in the radiation era. Matter era loops contribute to lower 
frequency part of the spectrum. 3 Since we want to get an estimate of the 

3 It is relatively easier to verify this in the case of small loops. If one limits the redshift 
integration in Eqs. (62) and (63) to matter domination and uses the corresponding ap- 
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spectrum we will neglect the complications arising from removing rare 
burst. In the radiation domination, z > z eq = \ATr — 5440, the Hubble 
function in Eq. (33), can be approximated as 

H(z) ~ yftte? = (58) 



The cosmological functions that appear in the stochastic background 
radiation formula can be approximated as 

/ x f°° dz' f°° dz' o , N 

VAZ) = 1 (1 + z')H(rJ) * I 7W7)~^ Z ~ ' (59) 

, . f z dz' f z "i dz' f z dz' , . 

^ z) = iw) = l W) + Lwr^ m 

^ 1= (H%W gig ^ (61) 
We first consider the small loop case, for which the expression for SBGW 
reduces to an integral over redshift given in Eqs. (51) and (52). Inserting 
the result in Eqs. (59-61) into Eq. (51) we get 




H z 



-1/3N 



oc ^, (62) 
V 



where we dropped a term with 1/ / dependence since it is small in large / 
limit, and the subscript R reminds us that this is the contribution from 
radiation era loops. The upper limit of the integration, z max , is the redshift 
at the time of the creation of the strings, which depends on the energy scale 
of the phase transition. The result in Eq. (62) is valid for 

z 1/2 f 

<C -ft- < -^ffk, for which the upper limit of the integral is set by the 

OLZ e q 

O-function. If > -^ffk, the integral does not depend on / and the 



proximate cosmological functions, it is found that £l(f) depends on the negative powers 
of /, which are negligible for large /. The same argument also applies to the large loop 
case. 
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frequency dependence of Q gWt R,(f) is given by the prefactor, which has / 1 ^ 
behavior. For kinks we get 



o K m nr Gfl f Zmax dz n I 1 



#0^ 



-1/3N 



oc^. (63) 
P 



Eqs. (62) and (63) show that for <C ^- the spectrum is constant and it 
scales with G/i/p. The amplitude does not depend on the parameter a, 
however the spectrum shifts to the right linearly in a. 
This result is in perfect agreement with Fig. 1. For the bottom curves 
^ = 2x 1CT 5 where as — = 2 x 10" 3 for the top curves, which have two 
orders of magnitude larger amplitude, exactly agreeing with the figure. 
Furthermore, the top curves (e = 1CT 4 ) are shifted to the right compared to 
the bottom curves (e = 1) by about 4-orders in / as predicted by our 
results above. 

Now we consider large loops in the radiation domination, for which the 
density n(l,t) is given in Eq. (53), where t is to be replaced with (p t (z)/Ho. 
Substituting the results in Eqs. (59-61) into Eq. (38) we get 



(lz 2 + (3S)l fzl 

l -A 

A(f) I dz j z du U - g (64) 

2e<7 J J (UZ + (35) 2 



where we define 



A(S) = ^^'f^, , (65) 

with 5 = and (3 = (a ~ 0.1 for large loop case) and the dummy 

integration variable u = I z. The upper limit of the z integral, z* will be set 
by requiring (3/z > 1/f, that is, z < ff3. If / < z max /f3 we have, 

n(/) = A(f) [ m dz [ § du u ~ g =a(/) /^d« r<fe u ~* 5 



2 . , „ s fi^ ofai / 1 1 



-^ (/) I 1 ul{((3 + /3*)§ (uz eq + 0S)i) ' (66) 
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If j < y~ — g ^E_ , we can split the integration range [1/ f, /3/z eq ] in the 

second integral into [1/ f,Sj3/z eq ] and [5f3/z eq , (3/z eq \ and neglect uz eq and 
(35 respectively in these two integrals. Combining all terms and keeping the 
lowest order in 5 we get, 



fW/) = A(f) 



2/3 18z eg 3 \ 330 c a 2 5-2 X r 



(6/3)1 ll(Sf3) 11/6 J pZeq i/*H$r*f3i 



* 3.2xl0-^/> 3 - 6 ^ 10 " 18 ^ (67) 

The calculation for the case of kink is very similar to cusp case, following 
the same steps we get 

flJU/) - 3.2 x 10- 4 ^, / > ™*™~™ Hz (68) 

which is identical to the cusp result. Eqs. (67) and (68) show that the 
distribution is flat for / > 3,6 — Hz and its amplitude scales with 
y/Gfi/p, which is in excellent agreement with Fig. 2. The flat value of the 
spectrum for the top curves (Gfx = 1CT 7 and p = 5 x 1CT 3 ) is 2.1 x 1CT 5 
and for the bottom curve (G/u, = 1CT 9 and p = 5 x 1CT 2 ) is 2.1 x 1CT 7 . 
These results are to be compared with the analytical results 2.0 x 10~ 5 and 
2.0 x 10~ 7 predicted by Eqs. (67) and (68). 

It is important to note that, in this paper we assume that the number of 
kinks, N, is order of one. This assumption enters in the estimation 
\X±\ ~ y-, and if there are iV kinks on strings, it needs to be replaced by 
|^±| ~ jf^- The replacement of / with l/N should also be done in the 
opening angle of the cone of the radiation, Eq. (37), which will result in a 
nontrivial dependence on N. However we can simply convert the resultant 
expression to the one we calculated in Eq. (52 ) by defining a = ct'N. Since 
we have shown that a has the effect of moving the spectrum horizontally, 
one effect of having iV kinks will be shifted spectrum compared to one kink 
spectrum. The other effect will be an overall scaling of the spectrum by 
l/N. 
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5 Parameter Space Constraints and Results 

In this section we discuss certain experimental bounds on SBGW. For the 
case of large loops the parameters are G[i and p, and for small loops the 
parameters are G/i, e and p. It is important to note that the nontrivial 
dependence on p follows from excluding rare bursts as described in Eqs. 
(56) and (57) (if rare events were included would simply scale with 

VpO 




Figure 3: Top-left: Accessible regions in the e — Gfi plane for p = 10~ 3 for 
small loops (loop sizes are determined by gravitational back-reaction). Top- 
right: Same as above for p = 1CT 2 . Bottom-left: Same as above for p = 1CT 1 . 
Bottom-right: Accessible regions in the p — Gfj, plane for the large long-lived 
loop models. The accessible regions are to the right of the corresponding 
curves. All models are within reach of LISA and advanced LIGO, and most 
are within the projected pulsar bound. 
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Accessible regions corresponding to different experiments and bounds are 
shown in Fig. 3. The shaded regions, from darkest to lightest, are: LIGO 
S4 [33] limit, LIGO S5 [34], LIGO H1H2 projected sensitivity 
(cross-correlating the data from the two LIGO interferometers at Hanford, 
WA (HI and H2)), and AdvLIGO H1H2 projected sensitivity. All 
projections assume 1 year of exposure and either LIGO design sensitivity or 
Advanced LIGO sensitivity tuned for binary neutron star inspiral search. 
The solid black curve corresponds to the BBN [35] bound, the dot-dashed 
curve to the pulsar bound[36], the +s to the projected pulsar sensitivity, 
the circles to the bound based on the CMB and matter spectra [37], the xs 
to the projected sensitivity of the LIGO burst [22] search, and the o-curve 
to the LISA projected sensitivity [38]. The BBN and CMB bounds are 
integral bounds, i.e. they are upper limits for the integral of over In/, 
therefore a model is excluded if it predicts an integral larger than the limit. 
On the other hand, the pulsar and LIGO bounds apply in specific frequency 
bands, thus a model is excluded if it has £l(f) larger than the limit (or 
projected sensitivity) for any / in the range of the pulsar or LIGO 
experiments. The range of the redshift integral in Eq. (38) must chosen 
properly for a given experiment. For BBN bound, the integration is 
performed for z > 5.5 x 10 9 . Similarly, for the bound based on the CMB 
and matter spectra, the integration is performed for z > 1100. First, we 
note that smaller values of p are more accessible, which follows from the 
fact that the loop density is inversely proportional to p. This makes cosmic 
superstrings more accessible than field theoretical strings. Second, we note 
that LIGO stochastic search constrains large Gp, small e part of the 
parameter space, whereas pulsar limit constrains large Gp and large e part 
of the parameter space. Similarly, the LIGO burst bound applies to large 
Gp and intermediate e part of the parameter space. Therefore large Gp 
part of the parameter space is covered by these three experiments. 
Furthermore since they also overlap for large Gp and intermediate e, in the 
case of detection, the two LIGO searches could potentially confirm each 
other. We also see that the BBN and CMB bounds are not very sensitive to 
e: the corresponding curves are rather vertical in e — Gp plane. This result 
is in perfect agreement our results (Eqs. (62) and (63)) that show 
fi(f) oc Gp/p, which does not depend on e. For the case of large loops, GW 
background is significantly larger than the small loop one, see Figs. 1 and 
2. Therefore more of the parameter space is accessible to the current and 
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proposed experiments, as depicted in the right bottom panel of Fig. 3. The 
strongest constraint is the pulsar bound, which rules out cosmic 
(super)string models with Gfx > 1CT 12 and p<8x 1CT 3 . This bound also 
rules out field theoretical strings (p = 1) with GpL > 2 x 1CT 9 . One can 
compare these results with the case where only cusps are included [25]. In 
that case cosmic (super)string models with Gfi > 10~ 12 and p<3x 1CT 3 
and field theoretical strings with Gji > 10~ 9 are ruled out. This result 
illustrates that kinks contribute to SBGW at the same order as cusps. 
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